Thermal Logic Gates: Computation with phonons 
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Logic gates are basic digital elements for computers. We build up thermal logic gates that can 
perform similar operations as their electronic counterparts. The thermal logic gates are based on 
the nonlinear lattices, which exhibit very intriguing phenomena due to their temperature depen- 
dent power spectra. We demonstrate that phonons, the heat carriers, can be also used to carry 
information and processed accordingly. The possibility of nanoscale experiment is discussed. 
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As two fundamental energy transport phenomena with 
similar importance in nature electric conduction and 
thermal conduction have never been treated equally. 
The invention of electronic transistor Q and other rel- 
evant devices that control electric flow has led to an 
impressive technological development that has changed 
many aspects of our daily life. However, similar de- 
vices in controlling heat flow are still not available even 
though many experimental attempts have been made 
to design such devices 0, 0, IE & 0|- Recent years has 
witnessed some important progress, for example, ther- 
mal rectifier has been proposed by using the property 
of temperature dependent power spectra in nonlinear 
lattices]!, S ES Ell, E3, H Q O , and molecular level 
thermal machines that pump heat from a cold to a hot 
reservoir has been suggested lfl 17 1. Merely few years af- 
ter the theoretical models of thermal rectifier, a nanoscale 
solid state thermal rectifier has been demonstrated exper- 
imentally by asymmetrically deposited nanotubes [18j |. 
Moreover, based on the new phenomenon: negative dif- 
ferential thermal resistance @, we have also built up a 
thermal transistor modelpjj which controls heat flow like 
a Field-Effect- Transistor (FET) does for the electric cur- 
rent. 

On the other hand, phonons and other vibrational ex- 
citons in low dimensional adsorbed dielectric and semi- 
conductor nano structures might also play an important 
role in information transmission (20j . 

In this Letter we present thermal logic gates that can 
do basic logic calculations. This may provide the pos- 
sibility that in the near future, phonons, which are tra- 
ditionally regarded as the heat carriers, can also carry 
information and can be processed accordingly, like elec- 
trons and photons, 'thermal computer', an alternative to 
the existing electronic computer, may also be possible. 

In an electric/chemical[2lj logic gate, power supplies 
that fix the voltages/chemical concentrations of some 
points are necessary. Similarly, in a thermal logic gate, 
power supplies that keep the temperatures of some points 
are necessary too. The temperatures of the power sup- 
plies are labeled as T + and TL, T + > T_. In any lin- 
ear electric circuit in which all the resistances are con- 



stants, when the voltage of one node is changed by a 
battery connected to it, the voltage of any other node 
may also change, but must be in the same way, and the 
latter change must NOT be greater than the former one, 
namely, 'super response' and 'negative response' axe not 
allowed. This is also true for thermal circuits. There- 
fore, nonlinear devices are necessary for any thermal logic 
gate. 

The nonlinear device that we use in this Letter is the 
thermal transistor, which is depicted in Fig[TJa)). Two 
weakly coupled nonlinear segments, D and S, are con- 
nected to power supplies with temperature T + (=0.2) and 
T_(=0.03), respectively. The control segment G is cou- 
pled to the 'input signal' with temperature Tq. Each 
segment is modeled by a one-dimensional (ID) Frenkel- 
Kontorova (FK) lattice [12] whose Hamiltonian reads: 



H 



FK 



x^f + U i {x i )]. (1) 



where k is the spring constant and Ui(x) = 1 — 
? cos 2ttx corresponds to the on-site potential. The 
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FK model describes a chain of harmonic oscillators sub- 
ject to an external sinusoidal potential. This is similar 
to the case by putting a p olymer chain or a nanowire on 
top of adsorbed sheet [2^. Since the momentum conser- 
vation is broken in such a lattice, heat conduction obeys 
Fourier's law 23|. For simplicity, we have set the masses 
of all particles be unit. The spring constant k and the 
on-site potential strength V in different terminals are dif- 
ferent. Terminal D is coupled to terminal S by a spring 
of constant hint, and terminal S is coupled to control 
terminal G by a spring of constant fc, n t G . ki n t is the 
most important parameter of the thermal transistor. It 
controls the magnitude and position of the 'negative dif- 
ferential thermal resistance' effect discussed below. More 
details can be found in Fig 3(a) in ref.[19(. 

In our computer simulation, power supplies and input 
signals are simulated by Langevin heat baths and we have 
checked that our results do not depend on the particular 
heat bath realization (e.g. Nose- Hoover heat baths). We 
integrate the differential equations of motion by using the 
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5th-order Runge-Kutta algorithm. The simulations are 
performed long enough to allow the system to reach a 
non-equilibrium stationary state and the local heat flux, 
Ji = k(xi(xi — Xi-x)), is constant in each segment. Tem- 
perature of a particle is defined as the average of twice 
of the kinetic energy, Ti — (x?). 

If the lattice in segment D is linear, then as tempera- 
ture difference in this segment is decreased by increasing 
Tq, one expects an decrease in J D . However, the FK lat- 
tice is nonlinear, one actually obtains an increased Jd as 
is demonstrated in FigQJb). In this figure, when Tq (Tq 
is always very close to Tq) changes from T ff(= 0.03) 
to T on (— 0.14), the temperature difference Td—Tq is de- 
creased, but the current Jd increases from 2 x 10~ 6 to 
5.9 x 10~ 4 which is about 300 times. This is the so- 
called negative differential thermal resistance (NDTR), 
an essential physical principle for thermal transistor [19j. 
The NDTR can be understood from studying the vibra- 
tional spectra of interface particles O and O' . One should 
recall that in our model, the heat is conducted by lat- 
tice vibration. Therefore, when one combines two seg- 
ments of different lattices, the overlap of the vibration 
spectra of the interface particles mainly determines the 
heat current. If the two spectra overlap, the heat can 
easily go through, otherwise, it is much more difficult. 
The FK lattice has a temperature dependent vibrational 
spectrum. At the high and low temperature limits, the 
power spectra concentrate on low (o;£E[0, 2-\/&]) and high 



the transistor shown in FigQJa) 



(lu e[vF, y/V + 4/c]) frequencies, respectively. More de- 
tails can be found in Ref. [3] ■ As is shown in the inset 
of FigQJb), at Tq = T ff, the spectra of O and O' do 
not overlap, thus the heat current at this point is very 
small, whereas at Tq = T on , the overlap is much better, 
thus although temperature difference at the interface is 
much smaller, the heat current is however much larger. 

NDTR provides the keys for thermal logic gates, e.g., 
'super response' and 'negative response', which are nec- 
essary for thermal signal repeater and NOT gate. See 
Figf2{a), in the central region, when temperature Tq 
is changed, temperature To changes even greater than 
Tq, which makes To be always closer to T on or T" // 
whichever is closer, than Tq is. This can be understood 
by considering the direction of heat flow in segment G 
when Tq is slightly different from T on /T ff. And see 
FiglUb). To> changes in the opposite way of Tq. This 
can be understood by considering change of heat flow Jd 
versus To- As Tq increases (in this case the thermal re- 
sistance between G and O is very small thus To is always 
close to Tq) thanks to the NDTR effect, heat flow Jd in- 
creases, which makes temperature difference Td-To' in- 
crease, namely To> decreases. As have been mentioned 
above, these 'super response' and 'negative response' are 
not possible in any linear circuit. In the following, we 
show how to build thermal gates by combining thermal 
transistors. Hereafter in this paper, symbols O, O' , D, S 
and G particularly indicate the corresponding nodes of 
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FIG. 1: (Color-online) (a) Configuration of the thermal tran- 
sistor, (b) Heat currents through three terminals D, S and 
G versus temperature Tq. Notice the NDTR effect that in a 
wide region both Js and Jd increase when temperature Ta 
is increased. At Ta—T on and T //, Jq is exactly zero. In- 
sets: power spectra of particles O and O' near 'off' and 'on' 
states. Power spectrum of O depends on temperature sensi- 
tively. It matches that of O' much better at 'on' state than at 
'off' state, thus makes much higher Jd at 'on' state although 
temperature difference between terminal D and particle O is 
much larger at 'off' state. 

In a digital electric circuit, two standard values of volt- 
age are used to indicate states '1' and '0'. Similarly we 
use two standard values of temperature T on and T a f f and 
hereafter we fix T o „=0.16 and T o //=0.03. 

The most fundamental logic gate is the signal repeater, 
which is a two-terminal device (one input and one out- 
put) whose function is to standardize an input signal. Its 
response function, T ou t P ut as a function of Ti nput , is: 



Toutput — off, 

Toutput — T on , 



if T, 



input 



if T, 



input 



<t c 

>T C . 



(2) 



Namely, when input signal is lower/higher than a critical 
value T c (T ff<T c <T on ) the output is exact T ff/T on . 
This is not a trivial device, without which small errors 
may be accumulated thus eventually leads to a wrong 
calculation. 
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FIG. 2: (Color-online)(a)r -T G versus T G . Notice the 'Super 
Response' that in the central region To-Tg increases when 
Tg increases, namely To increases even greater than Tg does, 
(b) Tqi versus Tg- Notice the 1 Negative Response" that Tqi 
decreases as Tg increases. 
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FIG. 3: (Color-online) (a) Structure of a thermal repeater, 
(b) Function of a thermal repeater that contains six thermal 
transistors. The outputs are better and better as the number 
of transistors is increased. The final output from the 6th 
transistor, is very close to that of an ideal repeater. 



This ideal repeater can be realized by several two- 
terminal devices whose response functions have two stable 
fixed points at T on and T Q ff. Therefore, if we connect 
them in series, namely plug the output of one device 
to the input of the next one, the final output is closer 
and closer to an ideal repeater. It is easily seen that 
in such a device, 'super response' i.e., T ou t P ut changes 
even greater than Ti nput does is necessary, otherwise the 
two fixed points at T on and T Q ff cannot be both stable. 
Such a device cannot be realized by any linear thermal 
circuit however can be simply realized by a thermal tran- 
sistor by using terminal G as input and node O as output 
(Fig|2^a)). Tq = T on and T a ff are the two stable fixed 
points and there still exists an unstable fixed point at 
Tq = T c , which separates the attraction basins of the 
two stable fixed points. 

As an example, we show a repeater consists of six ther- 
mal transistors in FigG^a). Its function is shown in Fig. 
[2Kb) . It is very close to an ideal repeater. 

A NOT gate reverses the input, namely if 
Tinput— T on /T ff then T ou tput—T f / /T on . This re- 
quires that when Ti nput increases T ou tput decreases, and 
vice versa, which corresponds to the 'negative response', 
which is, as mentioned above, again realized by a thermal 
transistor by inputting signal to G and then collecting 



output from O' , see Fig{2f b) . Notice that To< is always 
higher than T c (in fact even higher than T on ), thus will 
be always treated as 'on' by the next device. In order 
to solve this problem we use a 'temperature divider' 
(counterpart of a voltage divider which is shown in inset 
of FigH] (b)) whose output is a ratio of its input. By 
adjusting this ratio, we make its output be higher/lower 
than T c when the original input Tq is T ff/T on . Then 
after being standardized by a thermal repeater the final 
output realizes the function of a NOT gate. See FigQ] 
(a) and (b) for its structure and function. It is quite a 
surprising result because high temperature means more 
entropy or more randomness. Here we demonstrate that 
one can reduce the randomness of one part of the device 
by inputting more randomness in another parti 

An AND/OR gate is a three-terminal (two inputs and 
one output) device. If one input is 'off '/'on' then the out- 
put must be 'off'/'on' regardless the other input. Other- 
wise if one input is 'on'/'off' then the output follows the 
other input. By plugging two inputs (or better after stan- 
dardized by a repeater) to the same repeater, it is clear 
that when both the two inputs are 'on'/'off' the output 
must be 'on'/'off'. By simply changing some parameters 
of the repeaters, we can make the final output 'on'/'off' 
when the two inputs are different, therefore AND/OR 
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FIG. 4: (Color-online) (a) Structure of the thermal NOT gate. 
Through the G segment, signal is transferred to O point of 
the transistor. The output of the transistor (from O' point) is 
transferred to the temperature divider. Plug the output of the 
temperature divider into a repeater the final output is shown 
in (b): function of the thermal NOT gate. It is very close to 
an ideal NOT gate. Inset: Structure of a two-resistor voltage 
divider, the counterpart of a temperature divider, which sup- 
plies a voltage lower than that of the battery. Without load, 
the output of the voltage divider is: V ou t = VR,2/(Ri + B.2)- 
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FIG. 5: (Color-online) (a) Structure of thermal AND and 
OR gates. Two inputs (after standardized by repeaters) are 
transferred to a same repeater, (b) Function of thermal AND 
gate. We plot T outp ut versus T inp ut2 while T input i is fixed to 
T on . (c) Function of thermal OR gate, we plot T outpu t versus 



T, 



in P ut2 



while T inp uti is fixed to T D /j 



gate is realized. See FigJSJa) for their structures and (b) 
and (c) for their functions. 

In summary, we have presented the feasibility to realize 
thermal logic calculation based on the thermal transistor 
model. Although the thermal logical gates given here 
are only toy models, they may shed light on the study of 
molecular information technology [2J| and smart thermal 
materials. The study may also be helpful in understand- 
ing the complicated heat transport in biological systems. 

The FK model we have used here is a very popular 
model in condensed matter physics and nonlinear physics 
22] . It is usually used to describe energy transport in real 
solid. We therefore believe that our thermal logic gates 
can be realized in nanoscale systems experimentally in 
a foreseeable future, in particular, given the fact that 
the solid state thermal rectifier has been realized experi- 
mentally in 2006 [l8j], which is only a few years after the 
theoretical works. 

The work is supported in part by an Academic Re- 
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